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Abstract 

We consider self-similar (pseudo-steady) shock reflection at an oblique 
wall. There are three parameters: wall corner angle, Mach number, angle 
of incident shock. Ever since Ernst Mach discovered the irregular reflec- 
tion named after him, it has been an open problem to predict precisely 
for what parameters the reflection is regular. Three conflicting propos- 
als, the detachment, sonic and von Neumann criteria, have been studied 
extensively without a clear result. 

We demonstrate that the sonic criterion is not correct. We consider 
polytropic potential flow and prove that there is an open nonempty set of 
parameters that admit a global regular reflection with a reflected shock 
that is transonic. 

We also provide a clear physical reason: the flow type (sub- or super- 
sonic) is not decisive; instead the reflected shock type (weak or strong) 
determines whether structural perturbations decay towards the reflection 
point. 



1 Introduction 



1.1 The transition problem in shock reflection 

Reflection of an incident shock from a solid wedge is a classical problem of gas 
dynamics. It has been studied extensively by Ernst Mach [211 [18] a-^d John von 
Neumann [25], among others. 

Most commonly, reflection is studied in steady inviscid polytropic0 compressible 
flow, for example when shocks in a nozzle are reflected from the walls. The 
reflections can be classified roughly into regular and irregular reflections; see 
[1] or [2] Figure 1] for a more detailed discussion. In either type, an incident 
shock impinges on a solid surface. In regular reflection (RR), the incident shock 
reaches a reflection point on the surface, continuing as a reflected shock (see 
Figure [1] top left). 



^equation of state p = (7 — l)pe, e internal energy per mass, 7 G (1, 00) 
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Figure 1: Left top: local RR. Left bottom: subsonic case. Right: fixed V2\ each 
steady shock produces one 1)3 on the curve (shock polar, symmetric across V2; 
shock normal || V2^vz). For |t| < t*, three shocks satisfy t = /(w2, v^)'- strong- 
type (S), weak-type (W) and expansion (U; unphysical). W are transonic right 
of supersonic left. 
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Figure 2: Left: initial data (from dotted area of Figure [3] second right). Center: 
RR; we construct perturbations of the trivial case. Right: SMR 
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Figure 3: A Shockwave breaks into a supersonic RR at the lower corner; the 
reflected shock breaks in the upper corner to produce our kind of subsonic RR 
for some time. 
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In irregular reflections (IRR), incident and reflected shock are connected by a 
more or less complex interaction pattern which in turn connects to the solid 
surface by a third shock, called Mach stem. The most important irregular 
reflections are double, complex and single (see Figure [H right) Mach reflection 
(MR); additional types have been discussed [H Hi [H [T71 US] . 

Some incident shocks allow more than one type of reflection. Assuming unique- 
ness for the problem at hand, only one of them can be extended to a global 
solution: a solution in the entire domain, satisfying all boundary and far-field 
conditions. A long-standing open question is to find the exact criterion that 
determines whether the solution is RR. 

Among the criteria for appearance of RR that have been proposed (see [1] 
Section 1.5]), two are most important. The detachment criterion states that 
global RR appears generically whenever local RR is possible. 

A physical argument motivates the second criterion: for a straight wall, all local 
RR and MR are trivially global solutions. But some of them could be unstable 
under perturbations, for example a bump in the 3-sector wall (Figure [T] left 
bottom). If so, then information is transmitted from the bump to the reflection 
poinfl For weak wave^ that is possible if and only if the 3-sector is subsoni^ 
Hence the sonic criterion: global RR appears generically if there is a supersonic 
local RR, but not otherwise. (Each criterion can also be formulated in other, 
slightly different ways.) 



1.2 Weak- and strong-type 

The velocity U2 in the 2-sector in Figure [T] forms an angle r with the wall; the 
reflected shock must turn this velocity by r so that W3 is parallel to the wall, 
satisfying a slip boundary condition. 

Keep the 2-sector data fixed while rotating the reflected shock in the reflection 
point. This yields a one-parameter family of velocities W3, forming a curve called 
shock polar (see Figure [T] right). For admissible shocks, |t| cannot exceed r*, 
the critical angle, which is a function of the Mach number M2 and 7. 

Throughout this paper we focus on polytropic equations of state so that the 
admissible part of the shock polar is strictly convex. 

If the angle r between wall and V2 is bigger than r, , then local RR is theoretically 
impossible. If r = t, , there is exactly one reflected shock, called critical-type. 

■^this is known as information condition or information argument 

^hut sufficiently strong shock waves can travel upstream against a supersonic flow 

*the other sectors are always supersonic 
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For T < T,t however there are two, called weak-type^ and strong-type. We en- 
counter another major issue in reflection: which of these two should occur? [lO] 
have discussed this question for a different problem. 

We call shocks transonic if the downstream side is subsonic, supersonic if both 
sides are supersonic. The weak shock is transonic for r > Tg for some threshold 
Ts < T», supersonic otherwise; the strong-type shock is always transonic. In this 
paper we consider only transonic RR. 



1.3 Self-similar reflection problems 

Some variants of the reflection problem are self-similai^ flow rather than steady. 

In self-similar flow, density and velocity are functions of the similarity coordi- 
nates (^,77) = {x/t,y/t) rather than x,y. Patterns grow linearly in time, with 
t I corresponding to "zooming infinitely far away" whereas t t -f-cxo is like 
"zooming into the origin" or "scaling up" . Here inviscid models are easily jus- 
tified because any flow feature eventually grows beyond the length scale where 
dissipate or kinetic phenomena mattei0. Self-similar reflections occur naturally 
in many experiments (see Figure [31 [T31[2]). 

Wc consider three parameters (see Figure [5] left): Mi, the 1-sector Mach are 
defined number, a, clockwise angle from opposite wall to incident shock, and 
180° — 6', clockwise angle from opposite wall to reflection wall. The opposite wall 
passeH through ^ ~ V2- Mach number and velocity are defined for an observer 
traveling in the reflection point. 

For t I this yields the initial datsfl seen in Figure [2] left. Depending on either 
RR or MR appear. 

If we choose the opposite wall perpendicular to the reflected shock, then local 
RR extends to a global trivial RR (see Figure [2] center) . 

^The weak-type shock is relatively weaker than the strong-type shock, but their absolute 
strength can be arbitrarily small or large, so we prefer to use -type. 
^also called quasi-steady or pseudo-steady 

^unless these small-scale phenomena trigger large-scale effects like turbulence, boundary 
layer separation etc. 

*to satisfy a slip condition on the opposite wall 

^If the incident shock forms a right angle to the upstream wall, this problem is familiar 
[21[8]. Note that the nonvertical cases also arise from certain t < flows; in particular they 
can arise in simple experiments like Figure |3] 
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1.4 Transition 



The two transition criteria specify transition angles 9d (detachment) and 6s > 9d 
(sonic) depending on Li, a, 7. Global RR is predicted for larger 6 and IRR for 
smaller 9. Figure |4] compares the two criteria in the case of 7 = 7/5 polytropic 
potential flow. 



To quote [J: 

"For this reason it is almost impossible to distinguish experimen- 
tally between the sonic and detachment criteria." 

Experimental and numerical accuracy are affected by viscosity/heat conduc- 
tioij^. non-equilibrum effects, turbulence, surface roughness and other system- 
atic errors as well as noise. The interaction of physical or numerical boundary 
layers with RR causes spurious Much stems [29i Figure 7a] that make it look like 
MR, however boundaries can be avoided by reflection into an interior problem. 



Although the question has remained open, the sonic condition appears to have 
been favored by many researchers (including the author), at least for small Mi. 
As the recent survey [2] states, 

"[...] the [criterion] which best agrees with pseudo-steady shock 
tube experimental data [...] suggests that in pseudo-steady flows 
RR terminates when the flow behind the reflection point, R [...] 
becomes sonic in a frame of reference attached to R." 

^"Observations (e.g. 1271 p. 142f) agree with our analytical solutions, so inviscid models 
arc clearly suitable. Experiments |14) show that, although viscous/boundary layer effects can 
have a transient effect on the transition 9, for sufficiently large times the transition is close to 
the inviscid predictions 9^,9^. 
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It should be noted that these quotes refer only to the classical case a = 90° 
and hence 6 < 90°, i.e. vertical incident shock (see Figure myreffigrexperiment 
second left). Here we consider some cases with 9 > 90° and a < 90° because 
they can be solved by linearization around trivial RR, i.e. small-data techniques. 
The classical case requires a large-data approach as in [T^; this will be subject of 
future research. However, the nature of the question is the same in all cases: how 
is local RR affected by various kinds of perturbation. The classical perturbation 
occurs naturally in some experiments, but there is no other reason to favor it. 



1.5 Results 

We prove, for the self-similar reflection problem modeled with potential flow, 
that the sonic criterion is not universally correct. We use the following formu- 
latioS 

Generic local transonic RR cannot extend into structurally stable 
global RR. 

Instead, Theorem [5] shows: 

1. Trivial weak- type transonic RR is structurally stable. 

In particular, the parameter space has an open nonempty — hence generic — 
subset with extendable local RR. 

More importantly, we identify a physical reason for the failure of the sonic 
criterion. The information argument (see above) indeed goes a long way 
towards the correct answer. But interestingly, it is too restrictive in a subtle 
way: 

2. For weafc-type transonic reflections, downstream perturbations 
can reach the reflection point, but they decay to zero^^ in the 
process. 

This suggests that — although a proof is given only for particular parameters 
— the sonic criterion is incorrect for most, if not all, parameters, in particular 
including the classical case a = 90°. 

We demonstrate the principle for a particular variant of the reflection problem 
in potential flow. However, it will become clear during the course of the proof of 
Theorem [2] that gradient regularity near the reflection point is a local property 
of elliptic PDE and their boundary conditions in a domain corner. Hence the 

^^This version is a weak as possible, by considering "generic" instead of all, and by requiring 
structural stability. 

^^Here we mean decay in space, not in time. 
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same principle applies to other steady or self-similar variants. Moreover, the 
same regularity effect occurs in isentropic and non-isentropic Eulcr flow. 

A wider range of parameters, Euler flow and the validity of the detachment 
criterion will be discussed in separate articles. A third important criterion, the 
von Neumann criterior^^ does not apply at all in potential flowF^. 

|10] previously provided a rigorous construction of supersonic weak-type re- 
flections in a difl'erent problem. The techniques in this article are apparently 
sufficient to extend the construction to transonic cases. 

Considerations analogous to the sonic criterion have also been used in studying 
the transitions between different types of Mach reflection (see the survey in [2] 
for details). Our findings suggest modifications to these proposals as well, by 
replacing subsonic-ness with weak-type in some way. 



1.6 Other remarks 

Many articles have considered dynamic stabilitAp^. trying to show that at the 
linear level weak- type reflections are stable while strong- type arc not. However, 
numerical calculations [D] Figure 3] suggest that both types are dynamically 
stable. 

[15] has previously proposed a plausibility argument, based on pressure changes, 
for stability of weak-type transonic (and instability of strong-type) shocks. [5] 
show existence and structural stability of supersonic reflections from a wedge. 
[3] have constructed global supersonic RR for a = 90° and 6 « 90° as exact 
solutions of self-similar potential flow. [5] shows existence of global supersonic 
RR for a range of parameters that includes, in some cases, all 6 > 6^, proving 
that criteria more restrictive than sonic cannot be universally correct. 

It has been proposed that both RR and MR may occur for the same parame- 
ters in steady flow, with hysteresis effects when parameters arc changed (see e.g. 
[5]). In self-similar flow this would amount to non- uniqueness for an initial- value 
problem. Indeed, [HIT] has found a set of initial data for the 2d Euler equations 
(both isentropic and non-isentropic) that appears to have two solutions, one the- 
oretical, the other clearly different and observed in all numerical calculation^^ 
For isentropic Euler, a rigorous proof of a different non-unic^ucness example has 
recently been proposed [TO] . 

^''also referred to as mechanical equilibrum criterion in some contexts 

^■'Even in Euler flow it applies only for large Mi, for example Mi > 2.2... for 7 = 7/5. 

^^stability under perturbations to the initial data 

^^In addition it is shown that the Godunov scheme can converge to either solution, depend- 
ing on the grid. 
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However, both results depend strongly on vorticity; uniqueness for the potential 
flow Cauchy problem is still expected and hysteresis is unlikely except as a 
transient phenomenon. 



2 Self-similar potential flow 

Here we prove technical results which are not previously available in the litera- 
ture. 



2.1 Equations 

2d isentropic Eulcr flow is a PDE system for a density field p and velocity field 
V, consisting of the continuity equation 

Pi + V • {pv) = (1) 

and the momentum equations 

{pv)t + V ■ (pu ® w) + Vp = (2) 

The pressure p is a strictly increasing smooth function of p. The sound speed c 
is 



V X 



If we assume irrotationality 
then we may take 

for a scalar potential (j). Assuming smooth flow, the momentum equations yield 

p = ^-i(A-0,-i|V0p) (3) 
where A is a global constant and where 



The remaining continuity equation ([T]) is unsteady potential fli 



ow. 



For any t we may change from standard coordinates (t, x, y) to similarity 
coordinates (t,^,??) with ^ = (Ci'?) A flow is self-similar if p^v are 
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functions of ^, rj alone, without explicit dependence on t. In potential flow that 
corresponds to the ansatz 

a^, 2/) = tip{x/t,y/t). 

By differentiating the divergence form ([1]) of potential flow and using ^ and 
(HI, we obtain the form 

{c'l - (VV- - C?) : V^^ = 0. (5) 

Here A : B is the Frobenius product trlA^B), :^ w w = wilF (not 'uFw) 
and is accordingly the Hessian. In coordinates: 

It is sometimes more convenient to use the pseudo-potential 
which yields 

(c^/ ~ Vx') : V\ + - |VxP = 0. (6) 
We choose ^ = so that 

P^^-'{-X-\\yx?)- (7) 

([6|) is manifestly translation- invariant. Translation is nontrivial: in (t,x^y) 
coordinates it corresponds to a change of inertial frame 

V ^ V — w, ^ ^ x/t ^ — w, 

where w is the velocity of the new frame relative to the old one. Obviously the 

pseudo-velocity 

z := Vx = VV- 

does not change. 

Self-similar potential flow is mixed-type; the local type is determined by the 
coefficient matrix c^/ — Vx^ which is positive definite if and only if L < 1, 
where 

^ \z\ \v — x/t\ 
c c 

is called pseudo-Mach number; for L > 1 the equation is hyperbolic. 
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2.2 Shock conditions 



The weak solutions of potential flow are defined by the divergence-form conti- 
nuity equation ([T|). Its self-similar form is 

V • ipVx) + 2p = 0. 
The corresponding Rankinc-Hugoniot condition is 

Pn< = Pdz'^i (8) 

where it,d indicate the limits on the upstream and downstream side and z", z* 
are the normal and tangential component of z. As the equation is second-order, 
we must additionally require continuity of the potential: 

tA" = V''- (9) 
By taking a tangential derivative, we obtain 

4 = 4 =: (10) 

Observing that cr = ^ ■ n is the shock speed, we obtain the more familiar form 

PuVu - PdVd = o{pu - Pd), (11) 

vi = v'd=:v'. (12) 

Fix the unit shock normal n so that > which implies > as well. To 
avoid expansion shocks we must require the admissibility condition z" > z^, 
which is equivalent to 

v: > v^. (13) 

We chose the unit tangent t to be 90° counterclockwise from n. 

By p2p the tangential components of the velocity are continuous across the 
shock, so the velocity jump is normal. Assuming w" > (positive shock 
strength), we can express the shock normal as 

- - (14) 



Wu - Vd\ 



2.3 Shock polar 



In our problem the upstream regions are constant and determined. Let ip be 
the potential in the downstream region, tp^ the potential upstream (ditto for x, 
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p, ...). We substitute HI]) into ^ to obtain the shock condition 

g(V^,V,a (^-'(-X- ^|VxnVx-/VxO ' |^^! I = 0. (15) 

The shock polar (see Figure [T|) is the curve of Vd that we obtain when holding 
the shock in a fixed ^ and keeping the upstream constant while varying the 
normal. For a fixed 5, Vx^ is fixed and i(j = ^/'(^) = ^^(C) is fixed as well. 
Having eliminated the normal in (fT5)) . we see that the shock polar is the curve 
of solutions V = Vtp of g{v,tp,S^) = 0. Hence the vector 

g- = (— — ) 

is normal to the shock polar, by the implicit function theorem. Omitting a 
positive scalar factor, it is given by the explicit formula 

g, ~ (1 - {z2/cf)n - z\\ + c-'4)l (16) 

as we show in 

For transonic shocks, which are our focus, the downstream is elliptic, i.e. 1 > 
Ld = \zd\lc > z2/c. In this case the coefficient of n in is necessarily 
positive, so g^; ^ 0. 

In Figure [T] right the leftmost point of the polar is a pseudo-normal shock: 
z* = 0. In this case (7,7 points in the same direction as n, hence right. Therefore 
g^; is an inner normaJ^J to the admissible part of the shock polar. 

In local RR the reflected shock must yield V3 parallel to the wall. In Figure 
[1] right, Vd for the weak shock (base in origin, tip in W) forms a blunt angle 
with inner normals of the shock polar whereas Vd for the strong shock (tip in 
K) forms a sharp angle. For the critical angle there is a single shock which is 
a limit of the weak and strong sides, so the angle is right (see in Figure [T] 
right). This motivates the following definition: 

Definition 1. A shock is called weak-type (in a particular point ^ in self-similar 
coordinates) if 

9v -ZdKO, (17) 
strong-type if negative, critical-type if zero. 



The definition has three pleasant properties: it coincides with the standard 
definition in the case of strictly convex polars, it generalizes the definition of 
weak/strong- type to non-convex caseJ^. and finally the sign condition is pre- 
cisely what is needed for elliptic corner regularity. 

^^not necessarily unit 

^*In such cases, there may be three or more reflected shocks that yield V3 tangential to the 
wall. 
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2.4 Polytropic pressure 



Throughout the paper we consider only the standard polytropic pressure law: 

1 Po 

with 7 e (l,oo), where co,po are constants. With this choice, 

Po 

Theorem 1. Consider arbitrary Cu,Pu > and Mu G (l,oo) and set = 
{MuCu,0). For each /3 G (—90°, 90°) there is a steady shock with downstream 
unit normal n ~ (cos/3, sin/3). Its downstream state pd, Cd, Vd depends smoothly 
on fi. Let t he counterclockwise angle from Vu to Vd- We restrict |/3| < arccos ^g- 
so that the shock is admissible. 

Then the shock polar f3 Vd is smooth and strictly convex, with dpVd nowhere 
zero. 

There is an angle r* G (0°,90°) so that each t G (— t*,t*) is attained for two 
different (3. The one with smaller \vd\ yields a strong-type shock, the other one 
weak-type. For |r| = they are identical and critical-type. 

There is a Tg £ (0,r^) so that the weak-type shocks are supersonic for \t\ > t^, 
transonic for \t\ < Tg. The other types are always transonic. 

Proof. We refer to [TU], especially Proposition 2.5.1, which establishes existence 
and smooth dependence of admissible shocks. By [TOl (2.5.2)] dpv 7^ at all (3. 

As shown earlier, g.g in (|16|) is an inner normal to the shock polar everywhere. 
Multiply it with a positive factor to obtain q ^ n — At where 

_ v\l/vl + Mll/cd) 

1 - {M2? ■ 

A is decreasing in /3 < 0, because by [ini Proposition 2.5.1] > is increasing, 
iij^ > is increasing, > is decreasing, ?;* > is decreasing. Hence dp A < 0. 
Moreover 

dpq ^ An + {1 - {dpA))t, 
{t is counterclockwise from n), so 

qxdpq = l- dp A + A^>0. 

This implies that the upper half of the shock polar is strictly convex. By vertical 
symmetry and smoothness the entire polar is strictly convex. 
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The shock polar is compact when adding the "vanishing" shock Vd = Cv ■ More- 
over |r| < 90°, so there is a maximum € (0°,90°). By convexity there are 
exactly two points on the polar for |t| < r*, which are the intersections of the 
line of multiples of Vd with the polar. As is an inner normal, necessarily 
gv'Vd>0 for the point closer to the origin (strong- type), with opposite sign for 
the other (weak- type). 



g, . (1 - K/cdfM - iv'ni/< + cfv-d) = v^dUl MJ) - 

If Md > 1, then the right-hand side is negative, so the shock is weak-type. By 
Proposition 2.5.1, Md is strictly decreasing in |/3|, so there is a unique r, so that 
the weak- type shock is transonic for |t| > t^, supersonic for |t| < r,. □ 



3 Perturbations of weak trivial RR 
3.1 Coordinate transform 

We consider a trivial RR as in Figure [5] center or Figure [5] left. All lines and 
curves exclude endpoints by default. We use the following notation (see Figure[6] 
left) : Let B be the reflection wall, A the opposite wall, W the open convex cone 
enclosed by them. Let nA,nB be the outer (with respect to W) unit normals 
of A, 13. Let the origin the the corner between A, B. Let S the reflected shock, 
£,A — (Ca,0), ~ {^b^Vb) (note S^A ~ £,b) the points where it meets A resp. 
B. Let A, B be the segments of A, B from the corner (0, 0) to £_aXb] let be 
the triangle enclosed by A,B,S. 

The velocity in is zero in the chosen coordinates, so the velocity potential i/' 
is constant = ■0° in Q,. Let vi — (w/,0) be the 2-sector, -0^ the corresponding 
potential. 

In self-similar flow, is a uniformly elliptic region whereas the rest of W is 
uniformly hyperbolic. 

The shock is a free boundary. To linearize the problem, we first devise a trans- 
form from ^ = (^, 77) to fixed coordinates a = (cr, 

Given a function e C^(ri) fl C^(r2). Consider a ray starting in the origin and 
passing through (^a,C) G S. is strictly monotone along any such ray, so 
there is a unique point ^ with 

^'(0 = V'(a,c)- (18) 
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(cr, C) enis mapped to 

This coordinate transform allows to state our problem in a, fixed domain Q. By 
(fT8|) . ip mapped to ^ coordinates satisfies the first shock condition, ip = , 
automatically. Then (fT5|) can be used as the second shock condition. 

3.2 Linearization 

We regard our problem as an operator equation 



where F : X ^ Y, X,Y Banach spaces with X C C^iQ) fl C^{n). The map 
is the composition of three steps: first tp X is transformed from (a, C) to ^ 
coordinates, then mapped to the tuple 



Finally, pull back this tuple to (cr, C) coordinates. 

X will be specified later since we have to consider an entire scale of such spaces. 
F will be a nonlinear map from X to Y. We intend to apply the implicit 
function theorem. To this end we need to study the Frechet derivative F'{ip'^) 
of F with respect to ip a.t ip = ip'^ . 

The derivative is computed by considering a first variation ip' ^ X oi ip*^. We 
consider = + tip'{(T) and compute the derivative ^ of F{tp) and other 

expressions, evaluated at t = 0. This derivative will be written p', F{'ipy, 
etc. Obviously the usual calculus rules apply. 

The following calculations are quite similar to [101 Proposition 4.14.3]. The 
results are simplified by two facts: ip = "0" yields an identity (^, 77) = (cr, (^), and 
= V2V'° = 0. 

All derivatives are evaluated at ip = tp*^ ; we omit arguments where they are clear 
from the context. 



F{ij) = 




[Interior] 

[Slip condition at A] 
[Slip condition at B] 

[Shock condition] 



V'ip\B ■ ns, 




(19) 
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=0 

^(i:S)V+E§(vF)' 

=0 

=0 =0 

A; 



Frechct derivative of the interior equation: 

= (c^/ - (V^-x)^)' : VfV' + {c'l (V^-x)^) : (V|^)' 

=0 

= {c'l ~ (V^-x)^) : V|^' = {c'l C") : V|^' (20) 

The resulting right-hand side is a linear eUiptic operator without zeroth-order 
term, applied to ip'. The classical maximum principle shows that ip' cannot 
have a minimum in the interior. 

The wall conditions linearize to 

VV'' • n = 0. (21) 



For the shock condition, we consider (jl5p . First, hold ^, ijj fixed and very 
The variation of the normal expression (jl4p is 



vi — VaI' 



-1 



VI - Vfipl \vi ~ V,-V'| 



1 - 



V^-^' = ^^^f (22) 
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Moreover, 

Both combined, we use the shock relations pj = pXn/Xn ''^^'^ Xt = xi to compute 
= p((l - c-\l)n -Xt{\+ c-\,,)t) ■ (23) 



Now we hold fixed and vary ip. ip = ij/ = t/j^ [0, 0) + vfS^ on the shock, so 



we can use 



= («?)-V'. (24) 
Moreover, the variation of the "normal" vj — is zero here, so: 

(5)' = (^"'( - ^ + ^leT - ^iVj-^ - cT)(v^-^ - a - - a)' 



= (1,0) 



""^ ' =(1,0) 

=v,-iA-(5)'=o 



-pC ^XnV'' + (P/ - P)C ^XnV'' + X 

Vj 

= -P{-^ +C~'^XnW 
Xn 

It turns out that rj' docs not appear in the final form, so the details of the 
coordinate transform do not matter at all! 

Altogether, when varying Vip and ip at the same time, the shock relations 
linearize to 

(ff)'=5.-V|^'-p(^+c-2x„)V'. (25) 

Xn 
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3.3 Kernel 



Proposition 2. For any X C C'^{n) n C^i^L), 

dimkerF'(?A°) < 1. 
If = \, then it is spanned by a function -0' that satisfies 

^'(6?) ^ 0. 

Proof. Assume the kernel is nontrivial. Let be a nonzero element. 

Consider a positive local maximum (with respect to Q) of i/;' at S* U {^a}- A 
maximum at S requires Tp'^ = 0; for a maximum in this is already implied 
by the boundary condition t/;^ = on A, by regularity in the corner, since 
A and S meet at a right angle. The coefficients of — V-'n a-nd -0' in the 
linearization of the shock condition, have opposite sign. Therefore ip' > in the 
maximum point implies —tpn > which is incompatible with a local maximum 
(n, the downstream normal, is an inner normal for Q). By the same argument 
a local negative minimum is ruled out. 

This implies in particular that tj/ cannot be constant. 

ip' satisfies pO|. the linearization of the interior PDE, so by the classical strong 
maximum principle ip' cannot have a local extremum in unless it is constant. 
By the Hopf lemma, the wall boundary condition ([21]) does not allow a local 
extremum at A or B unless "0' is constant. 

Assume ip' has a global maximum in (wall- wall corner). Let 3^(0) be the ball 
with radius e centered in and abbreviate U := B^{0) Ci fl, I := dB^{0) D il. 
For sufficiently small e > 0, I C flU AU B, so as shown above ip' cannot attain 
a maximum on /. Therefore tJj'^O) > maxjip' . 

tp :— ^' — ip'{0) + (5^ is a supersolution for S > 0: 

(I - c-2f ) : = (/ - c--'e) ■■ V^-' = 

by linearity, = on ^ and ipn = iS^)n > orF^ B. Therefore ip does not 
attain extrema in U. For 6 > sufficiently small, 

max0 = max0' - ^p'{0) + iJ^ > 0, 

while "0(0) — '0'(O) = 0, so the minimum of over U is attained in 0. Therefore 
0^(0) < 0, hence '0^(0) < —S < 0. But the boimdary conditions = on 

i^note e > 90° for trivial RR 
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A,B combine to V'0'(O) = — contradiction. Hence ip' cannot have a global 
maximum in 0; minima are ruled out analogously. 

Since ijj' is nonzero, it must have a positive maximum or negative minimum 
somewhere. As we have shown that is not possible except in S,r- 

For any two elements of the kernel, a suitable linear combination is zero in ^n, 
hence zero everywhere. Thus the kernel cannot have dimension higher than 
1. □ 



3.4 Type and Predholm index 

Proposition 3. Consider the eigenvalues of the operator pencil for F'{xp'^) in 
the reflection corner S,b (see Section \4.S^ . There is an eigenvalue Aq = ao + i(3o 
of multiplicity 1 with least nonnegative f3o, and 

{G (0, 1), if the shock is strong-type in ^b, 
= 1, for critical-type, (26) 
> 1, for weak-type. 

Proof. The operator F'^ipo)- with coefficients frozen in ^b, consists of the in- 
terior operator (/ — c^^^) : V^ip' and the boundary operators VV'' • n-B and 
gff ■ V'0'. We choose a linear coordinate transform so that the interior operator 
is mapped into Aijj' . This transform is a dilation in the B direction. 

Consider polar coordinates (r,(j)) centered in ^b- Let r2 = B, Ti ^ S, then 
the boundary operators take the form (PT|) with 72 = 90° (Neumann) and (see 
Figure [5]) 

{G (90°, 02 - 01 + 90°), for strong-type, 
= 02 - 01 + 90°, for critical- type, (27) 

G (02 - 01 + 90°, 180°), for weak-type. 

To see this, note that Vx° = V'0° — £, ^ \\ B on B. For a weak-type shock 
(Definition [1]) , g^j ■ Vx*' < 0, so ub x ga > 0. This property is preserved under 
dilation along B, so 71 > 02 — 0i + 90° (see Figure [5] left). The other types are 
analogous. 

Now (|32|) immediately implies p6|. □ 

Proposition 4. Consider the eigenvalues X = a + if] of the operator pencil of 
F' {ip^) in the A, B and A, S corner. The eigenvalue with least nonnegative (3 is 
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Figure 5: Corner limits of the top-order parts {gki,gk2) ■ V?/;' of the boundary 
operators. Left: weak-type shock. Right: strong-type shock: solutions need not 
be in the corner. A critical-type shock has (511, (712) exactly perpendicular 
to B. 



Ao = 0. The eigenvalue with next lowest nonnegative (3 is [3i = 1/(1 — 0/180°) > 
1 in the A, B corner and Pi ~ 2 in the A, S corner; their multiplicity is 1. 



Proof. In the A, B corner the interior operator is A, with Neumann boundary 
operators 9„ (so 71 = 72 = 90° in the notation of Section [42|l . so the calculation 
is straightforward. Take Fi := B, F2 := A. (pi = 6. 4)2 = tt, then by [32] 

/3o = 0, /3i = ^. 



In the A, S corner the slip condition on A yields Xv — 0; ^o the interior operator 
is (1 — x|)t^C5 + driri which becomes A by dilation. Moreover by ((25|) both 
boundary operators are which are not changed by dilation. Hence (5i — 2 
by □ 



Proposition 5. Consider a weak-type trivial reflection. Let s G (2,3). For 
e > sufficiently small, the Fredholm index of F'{ip'^) as a map from Xf^^ to 
Y-f^i^ (defined in Section \4-.2^ is 1. 



Proof. Consider the operator A — / on the triangle fio with Neumann boundary 
operators d„ on A,B,S. [101 Theorem 1.4] yields s £ (2,3) and e G (0, s - 2) 
so that the operator is a linear isomorphism on Xf+3g onto Yf_^_^i^. The space 

-X^i_l_3e defined in the present paper corresponds to ffi ^ ^'^^ in his notation, 
except that his weights are with respect to 9£2, not S. But d^l — E consists 
of line segments, so classical potential thcorjo [HI Lemma 6.27] extends his 
result to our case. 

■^''moreover s determines only regularity away from the corners, so it can be improved to 
any s > 2 
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an- 
vil' 
= 



n 

V^i^' = 

i' =0 



(A,S„) 



F'iii") 



Figure 6: Left: notation and linearized operator. Right: each point represents 
an operator T(*) : ^ Y| (t horizontal axis, P vertical axis). The operators 
are Fredholm except on the dashed curves, where some corner operator pencils 
has an eigenvalue. Between the curve the index is constant. Across the curves 
the index jumps by the eigenvalue multiplicity. 



Now we choose a family of operators t T^*) so that T^^-' ~ (A — /, 9„, 9„, (?„) 
and T^^^f = F'{iIj^) (see Figure |6]). The family is chosen continuous in t with 
respect to all operator norms we consider, which is easily achieved by choosing 
a continuous family of coefficients for interior and boundary operators. We 
choose the family so that /3o,/3i in the A,B and A, S corners arc constant in 
t (Proposition ID). If pf^ (j 0,1) are the two lowest nonnegativc imaginary 
parts of eigenvalues of the operator pencils in the corner, then t i-^ P^*^ are 
continuous as well. P^^^ = whereas Proposition [3] shows that /3q^^ > 1. By 
choosing suitable coefficient families in the reflection corner we can make Pq^ 
strictly increasing in t. Moreover p[°\ p[^'^ > 1, so we can achieve /jj*' > 1. 

m Theorem 6.3] yicldsE3 that T^*) : Xl_^^^ Yf^^ is a Fredholm operator if 
1 + e is not the imaginary part of an operator pencil eigenvalue in any corner. 

Choose e G (0,1) sufhciently small (not larger than above) so that 1 + 3e < 
P[*\l/{1 - e/l80°) for all t (see Figure H] right). Let e (0,1) be such that 

= 1 + 2e. Then t G [0,t,] ^ T^*) : F^g, and t £ 

T^*) : ^i+e ^''"6 both continuous families of Fredholm operators. By 

Their weighted Holder spaces are homogeneous; for our inhomogcncous spaces, ffilli (set 
of polynomials of degree < 1) is added for each corner, which is only a finite-dimensional 
change. 
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Fredholm theory the index of each family is constant. T*^") : — > ^i+e is an 
isomorphism, as shown above, i.e. has index 0. 

The interval [1+e, l+3e] contains only one eigenvalue of a corner operator pencil 
of T**, namely /Sj**"* = l + 2e; its multiplicity is 1. Hence [531 Theorem 6.4] shows 
that dim{Z/X+) = 1 in Proposition [SI where we choose X_ = Xf_^_^^, = 
Ff*+3„ X+ = Y+ = and A± = T^**). Therefore indyl_ - indyl+ = 1, 

so T(i) : ^ Yjt^^ has index 1. □ 

Remark 6. The proof requires /?o > 1, which is not satisfied for critical- or 
strong-type shocks in the reflection corner. This is the crucial difference to 
weak-type shocks. Note that the value of (3o is a purely local property; the 
chosen far- field perturbation is not significant. 



3.5 Perturbation 

Theorem 2. Consider a weak-type trivial transonic RR, with parameters po = 
{Ml, 6, a). There is a ball U of radius r > around po so that there is another 
global weak-type RR for any p £ U. 

Proof. Since the downstream state of a shock depends smoothly on the shock 
normal, location and upstream state, the shock polar varies smoothly with M„. 
Necessarily M2 > 1, so the incident shock is weak- type like the reflected shock. 
Therefore, sufficiently small perturbations p yield a new local RR which is close 
to the old one. In particular the perturbation of the reflection point is small 
and the reflected shock is still weak-type. 

By Proposition[5l the Fredholm index of F'{ip^) is 1. By Proposition[2]the kernel 
has dimension 1, so the codimension of the range is 0. Therefore we can apply the 
implicit function theorem, with a single real free parameter. By Proposition [2] 
we can use ipiS.s) as free parameter, which corresponds to changing the reflection 
point. 

Therefore we obtain a new elliptic region for sufficiently small perturbations of 
the reflection point, while satisfying both shock conditions. After extending the 
solution to the entire domain by adding incident shock and hyperbolic regions, 
we have obtained a global transonic weak-type RR. □ 



21 



4 Corner domains 



Here we adapt some literature results to our case. For details, see [12], [23] 
[22, and m. 



4.1 Weighted Holder spaces 

Consider a bounded open simply connected Lipschitz domain il C IR^. Let 
Ffc (fc = 1, . . . , m) be pairwise disjoint line segments with (excluded) endpoints 
yk-i,Vk (set yo ■.= y^^To := F„, for simplicity). Set E := {yi, . . . , j/^}. As- 
sume dVt = yy^^i Ffc. Let Fi, . . . , Tm pass around clockwise, so that lies 
counterclockwise from F^ to F^+i near each corner yk- 

Definition 7. Let /? e IR, s G (0, oo) - Z. Abbreviate fir := - -Br(S) where 
Br is the r-ncighbourhood. For u £ C^(r2,I]) we define weighted Holder norms 

ll"llcs(n,E) :=limsupr''-^||M||<.,(f^ ). (28) 

Then C^(r2, S) is the set of u with finite norm. The definitions for Fj in place 
of O are analogous. 

Non-integer (3 corresponds to the lowest exponent of r^ behaviour allowed in a 
corner; note that C|(n,S) C C'^(n). 



4.2 Operator pencils 

Consider the operator of a linear second-order elliptic boundary value problem: 

L(x)u :— y aijix)— — - — u+y bAx)— — u-\-c{x)u in 17, (29) 
axidxj ^-^ axi 

2 Q 

Bk{x)u:^'^gkt{x) — u + hk{x)u on Ffc, fc = 1, . . . , m. (30) 

i—i 

We assume that the coefficients aij,bi,c are smooth on il and gki,hk smooth 
on Ffc. We write B = {Bi, . . . , Bm)- Let s G (2, cx)) — Z; we use the convenient 
abbreviations 

m 

:= c^pin, E), c|:22(f], e) x [] c|il(r., e). 

fc=i 
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(L, B) : Xp -^Ypisa, continuous linear operator. 

Whenever L{yi) is elliptic, we can find a linear invertible coordinate transforma- 
tion so that the leading-order part of L{yi) transforms to the Laplace operator 
A. In this new frame we consider polar coordinates (r, (/)) centered in j/i. Let 
(j)i,4'2 correspond to Fi, we normalize 0i S [0,360°) and (?ii2 € [(j)i,(j)i+360°). 
The coordinate transformation from {x,y) to (t^(f)) with t = logr is conformal, 
hence preserves the Laplace operator, mapping the cone {(r,(l)) : r > 0, (p G 
(•/•ij </'2)} to an infinite strip IR x 02)- The leading-order parts of Bk{yi) are 

du du . 

_eos7. + ^sm7. (31) 

Here 7^ is the counterclockwise angle from to the coefficient vector {gui, 5fc2) 
on the corresponding boundary (see FigureE]). We normalize 71 G [0, 180°) and 
72 e (71 -180°, 71]. 



Apply the Fourier-Mellin transform in t to the homogeneous corner equation 

-l^(^t.,^)U = {-idtf + {-id^f = 
to obtain the operator pencil equation 

{—id^Yu + = 

where A = a + ifi are the eigenvalues. The eigenfunctions yield well-known 
harmonic functions 

u(i, 0) = exp(/3t) sin(/30 - 5) = r^ sin(/30 - 5). 

Imposing homogeneous boundary conditions restricts this family to 

u{r, 0) = sin{l3i{4> - (/)i) - 71) 

with 

02 - 01 02 - 01 

The multiplicity of each eigenfunction is 1, except when (3 = where it is 2 (for 
example in the case of two Neumann conditions there is another eigenfunction 
u = t = logr). 



4.3 Predholm index jump 



Proposition 8. Consider Banach spaces X+ C X„ and 1+ C Y- and Fredholm 
operators A± : X± '^Y±. Let Z -.^ {u X- : Au gY+] . If 

d :~ dim(Z/X+) < cxo, 
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then 



ind A_ — ind = d. 



Proof. ran^-|_ C ran^_, so kci A*_ — (ranyl^)-'- C (ranA+)^ = kcrA!j_; 
both spaces arc finite-dimensional by Fredholmness of A±. Choose a basis 
ipi, ■ ■ ■ ,ipr for ker C so that ipm+i , ■ ■ ■ jtpr form a basis for ker A*_ . Choose 
wi, . . . ,Wr € y+ biorthogonal to ipii ■ ■ ■ li'r- Then wi, . . . , Wm S (ker A*_)-^ = 
ranj4_ by choice of m, so we can find ui, . . . , Um G X- with Auj ~ Wj G y+, 
which also means ui, . . . , Um £ Z hy definition of Z . 

Claim: ui, . . . ,Mm,keryl_ are independent modulo X-^-. If not. we could find 
nontrivial coefficients ai, . . . , am as well as k £ kerA_, x e so that 

m 

QfiMi = k + X. 

j=l 

Then 

^ QfiMi = A/c + ^.T = Ax, 
i=l 

SO 

m m 

= V'jlX! "»'^') = V'il^X!"'"*) ^ = (j = 1, . . . ,n) 

i=l i=l 

since ipj G kerA'j_ = (ranA_|-)-'- and Ax E ranA_|-. The coefficients are trivial 
— contradiction. 

Assume we can add a u,n+i G ^ so that mi, . . . , Um+i, ker A_ are still indepen- 
dent modulo The system 

m+l 
i=l 

is underdetermined, so we can find a nontrivial solution ai, . . . ,am+i- Com- 
bined with the same result for j = m + 1, . . . , r (trivial) and with 

m+l 

'ijj{A ^ aiUi) =0 for t/i G kcr AI C ker A!j_, 

i=l 

we obtain 

m+l 

V'(A ^ Q!iMi) = for aU -ip G kcr 

i=l 

i.e. 

m-l-l 

A QfiUi G kcr(A!J.)^ = ranA_|_, 
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but then 

m+l 

A{Y^ aiu, -d) = 
1=1 

for some d £ so wi, . . . , Um+i, ker A- are dependent modulo — contra- 
dictfon. 

Hence ui, . . . ,Um form the basis of a complement of ker>l_ in Z modulo X^, 
dim(Z/X+) = dim(ker + to. 

Finally, 

ind A- — ind A^ 

= (dimker^_ — dimkerAl) — (dimker^+ — dimkcr^^) 

= (dimker^_ — dimker^_|-) + (dimkcryl!i_ — dimkcrAl) 

= dim(ker A_/ dimker74+) + to = dim(kerA_/X+) + to 

= dim(kcr^_/X+) + dim(Z/X+) - dim(ker = dim{Z/X+). 

□ 
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